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Assuming that the Hagedorn fluid composed of known particles and resonances with masses 
771 < 2 GeV obeys the first-order theory (Eckart) of relativistic fluid, we discuss the transport 
properties of QCD confined phase. Based on the relativistic kinetic theory formulated under the 
relaxation time approximation, expressions for bulk and shear viscosity in thermal medium of hadron 
resonances are derived. The relaxation time in the Hagedorn dynamical fluid exclusively takes into 
f^ , account the decay and eventually van der Waals processes. We comment on the in-medium thermal 

^\1 ■ effects on bulk and shear viscosity and averaged relaxation time with and without the excluded- 

volume approach. As an application of these results, we suggest the dynamics of heavy-ion collisions, 
non-equilibrium thermodynamics and the cosmological models, which require thermo- and hydro- 
dynamics equations of state. 






f-- ' PACS numbers: 66.20.-d, 25.75.Nq, 51.20.-|-d, 03.70.-|-k 

(N ■ 






o 



a 



I. INTRODUCTION 



Since the discovery of new matter at the relativistic heavy-ion collider (RHIC) [l|, an immense number of experi- 
mental and theoretical works have been invested to explore the properties of non-confined strongly coupled matter, the 
quark-gluon plasma (QGP). Characterizing QCD matter through transport and collective behavior and the equation 
Cn ■ of state is the ultimate goal. Also, the confined QCD matter has been subject of various papers [2t3- 

Shear viscosity r] characterizes the elliptic flow and is directly proportional to energy density e and inversely to 
scattering cross section a and simultaneously reflects how particles interact and collectively move in many particle 



> 

(«^ ' medium. Strongly interacting matter, like Hagedorn gas [6| is conjectured to have smaller r] than the weakly interacting 

^■f^ • one. Using perturbative and non-perturbative methods, t], (and normalized to entropy s) has been estimated for non- 
confined and confined QCD matter, respectively. 
At temperature T » Tc, bulk viscosity ^ has been studied using perturbative QCD Q and found to be negligible 

p^ ' comparing to r/ and regarding to the collective evolution of the many body system, as QCD turns to be conformal 
invariance. Recent lattice QCD simulations [7| show that bulk viscosity ^ is not negligible near Tc. Its rapid increase 
at Tc is apparently associated with a fast growth of the trace anomaly, (e — 3p)/T'^, of the energy-momentum tensor 
T'"^ Q. Below Tc, various hadron scales likely provide the conformal invariance with bad symmetries (QCD conformal 

S^ ' anomalous). It is therefore natural to expect that ^ at this energy scale is not negligible |10| . 

V^ I There have been many attempts to compute ^ and 77 in the confined QCD matter using effective models for the 
hadron resonances (a. Hll - llq . In this letter, we discuss the bulk and shear viscosity of the Hagedorn dynamical fluid 
using the relativistic kinetic theory in the relaxation time approximation 12tJ| and the explicit implementation of the 
hadronic mass spectrum and excludcd-volume approach. 

Disregarding all interactions but decay and repulsion, the thermal change of the relaxation time in the Hagedorn 
fluid has been studied, which complete a set of thermo [131 ^^"^ hydrodynamic equations of state needed to characterize 
the evolution equation in early universe, for instance. 

II. BULK AND SHEAR VISCOSITY 

The relativistic kinetic theory gives the transport equations for classic and colored particles in a non-Abelian 
external field. In this letter, we discuss the transport coefficients for confined QCD matter composed of known particles 
and resonances. We apply the relaxation time approximation of Boltzmann equation and take into consideration the 
hadronic density states in confined QCD. The transport equations describe the evolution of the phase space distribution 
function of the particles of interest. The transport properties are defined as the coefficients of spatial component of the 
difference between the energy-momentum tensors out of and at equilibrium corresponding to the Lagrangian density. 
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We model the Hagedorn fluid of QCD confined phases (fermion and boson resonances) as a noninteracting gas of 
hadron resonances. To do so wc sum over all Fermi and Bose resonances. The main motivation of doing this is 
that it refers to all relevant degrees of freedom of the confined strongly interacting matter. It implicitly includes the 
interactions that likely result in resonance formation [l8|. Natural units c = h = k = 1 are applied here. The Hagedorn 
mass spectrum p{m) implies growth of the hadron mass spectrum with increasing the resonance masses. 

p{m) = c (jtiq + 7Ti^) exp{in/TH), (1) 

with fc = -5, c = 0.5 GeV3/2, mo = 0.5 GeV and Th = 0.195 GeV. 

This model provides a quite satisfactory description of particle production in heavy- ion collisions [13, lU, [23] ■ The 
repulsive interactions likely soften the T-dependence of the thermodynamic quantities. The excluded-volume approach 
[l9[ is used to implement the effects of repulsive interactions (van der Waals) by assuming the energy normalized by 
43 equals the excluded-volume and the intensive quantity Tpt in the point-type particle approach (and the other 
thermodynamic quantities |20|) have to be corrected as follows. 

Tpt (^^^^ ^ M^l\ ^ (2) 



^ 4B \ ^ 4B / 

where B^/^ = 0.34 GeV stands for the MIT bag constant. 

In spherical polar coordinates, the energy-momentum tensor of a single particle with p- and T-independent mass 
m is defined as 

Tf = 7-^p(m) p^dp n{p,T), (3) 

27r^ J e 

where p^ = {e,p) is momentum four- vector and g is degeneracy factor of the hadron resonances. The single particle 
energy is given by the dispersion relation e = (p^ + ttt,^)^". 

With the above assumptions on Hagedorn viscous fluid, the overall energy-momentum tensor can be calculated as 
a sum over energy-momentum tensors T^'^ of all hadrons resonances ^ , 

rptliy ^ ^X^" (4) 

i 

In momentum phase space and assuming that the system is in a state with vanishing chemical potential and near 
equilibrium, the distribution function n{p,T) reads 

nip,T) = — ^^^ , (5) 

exp (^^) ± 1 

where ± stands for fermion and boson statistics, respectively. The local flow velocity u is compatible with the Eckart 
fluid [2j], implying that T^'^u^u^ = e. It is obvious that n{p,T) satisfles the kinetic theory [^ and second law of 
thermodynamics. The solution of kinetic equation is obtainable by deviating the distribution function from its local 
equilibrium. 

The deviation of energy-momentum tensor from its local equilibrium is corresponding to the difference between 
the distribution function near and at equilibrium, Sn = n — tiq. The latter can be determined by relaxation time 
approximation with vanishing external and self-consistent forces [3, UM 

Sn{p,T) = -T-r^df^72Q{p,T) (6) 

p ■ u 



Then the difference between near and equilibrium energy-momentum tensor reads 

9 

"27r2 



Srr - -T^Pi^) / P'dp^-T p"a„ no(p, T). (7) 



^ It reflects the algebraic properties, here the addition, of the energy-momentum tensor 



Using the symmetric projection tensor hajs [25| . the components of the derivative da can be spht to paraUel and 
orthogonal to u^. h^f^ generates a 3-matric and projects each point into the instantaneous rest space of the fluid. 

d^ = Du^' + V^, (8) 

where D = u^da = {dt, 0) gives the temporal derivative and V^ = 9^ — Ufj.D ~ (0, di) give the spacial derivative Q. 
Such an splitting has to guarantee the conservation of equilibrium energy-momentum tensor; dpT^'^ = and fulfill 
the laws of thermodynamics at equilibrium [3|, l25|. In ref [21, the non-equilibrium n{p,T) has been decomposed using 
the relaxation time approach n ~ nQ + rni + • • • . Alternatively, as n{p, T) embeds the l**-rank tensor m, 5r|"' can be 
decomposed into u [3| in order to deduce its spatial components. 

dku' = i Iduu' + 5nt' - '^SkAu^ + ^5kAu' = ]^Wm + ^5kAu\ (9) 



Applying the equation of hydrodynamics, then the deviation from equilibrium 

p'dp, (10) 



•^^r = ^^M'^ / ^^'^0(1 + no) 



p-UC^VaU +p — 7^^ ] +-:r-T\/aUi3 
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can be re- written as 



^T^ J Jo 
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rdp, (11) 



where c^ = dp/de is the speed of sound in this viscous fluid. The bulk ^ and shear 77 viscosity can be deduced from 
Eq. [11] by comparing it with 

STl' = ~^S,jdku''^7jW,j. (12) 

To find shear viscosity rj^ we put i 7^ j in Eqs. ([9]) and (|12p . To find bulk viscosity ^, wc substitute i with j and Tj^'' 
with 3P. The subscript (as in the distribution function n) refers to the equilibrium state. Although we keep the 
gradients of velocity, we put u = in the final expressions. The intensive quantities 77 and ^ of Hagedorn fluid ^ in the 
comoving frame read 

e(T) = ^l^p(m,)^°°no(l+no)(^c^£,?-ip2J p^dp (13) 

V{T) = ^LY^p{,n,)J^ no{l + no)^p'dp. (14) 

The T-dependence of dimcnsionlcss ratios (,T^/t and rfT /t is depicted in Fig. [1] With increasing T, bulk and shear 
viscosity increase, significantly. We note that ^ seems to be about one order of magnitude larger than r]. Left panel 
of Fig. El illustrates such a comparison. At low T, rj starts with larger values than ^'s. But with increasing T, ^ gets 
larger 12'^ . 

The ratio ^/ry has been related to the speed of sound Cg in gas of massless pions. Apparently there are essential 
differences between this system and the one of Hagedorn fluid. According to [23 , the ratio of ^/rj m. N = 2* plasma 
is conjectured to remain finite across the second-order phase transition. This behavior seems to be illustrated in 
Fig. [TJ In the Hagedorn fluid, the system is assumed to be drifted away from equilibrium and it should relax after 
a characteristic time r. Should we implement a phase transition in the Hagedorn fluid, then r oc ^^, where z is the 
critical exponents, likely diverges near Tc- 

III. RELAXATION TIME 

The relaxation time depends on the relative cross section as 

"(^) = nj[T){v{T)a{T)y ^^'^ 



As wc assume a vanishing chemical potential, the heat conductivity vanishes as well. 
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Fig. 1: Bulk (left panel) and shear (right panel) viscosity coefficients of Hagedorn fluid are depicted as function of the heat 
bath temperature. The effects of excluded- volume approach are illustrated (dashed lines). 



where v(T) and nf(T) is the relative velocity of two particles in case of binary collision and the density of each of the 
two species, respectively. The thcrmal-avcragcd transport rate or cross section is {v(T)a{T)). The transport equation 
of single-particle distribution function in the momentum space, n{r,p,t) [28|, 
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First line in r.h.s. of Eq. ([TS]) gives the Boltzmann collision term. The second line adds the Uehling-Uhlenbeck factors. 
The third line accounts for the Pauli-blocking of the final states. The total derivative of n is given by the collision 
integral. To solve Eq. (|T6|) . several gradients must be take into account. Real and imaginary parts of the G-Matrix [30| 
are taken to describe the potential of nuclear interactions and the cross section of the binary interaction, respectively. 
In-medium effects in final (Pauli blocking) and also in intermediate states have to be taken into consideration. In 
Ref. [26| . the kinetic Boltzmann-Uehling-Uhlenbeck equations for pure nucleon system have been analyzed and the 
relaxation time in non-relativistic approximation has been deduced as 
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where n is the baryon density and tlq is the nuclear saturation density np « 0.145fm~'^. 

When fitting the decay widths Ti of i hadron resonances, then the decay relaxation times ti in GcV^^ read ti = 
T~^ = (O.lSlrrii — 0.058)~ [13l. [20. |31| . As the resonance mass m is conjectured to remain constant in thermal and 
dense medium, this linear fit apparently implies that r remains unchanged as well. 

In the Hagedorn fluid, where the inter-particle collisions as in Eq. (|17p are minimized, we are left with specific 
processes to estimate r (decay and repulsion for instance). Formation from free space vacuum and decay to stable 
resonances; Pi+ P2 ^ P3 [32| are examples. The constrains on this process have been discussed in |32| . 

In rest of frame of the particle P3 boosting from the laboratory frame, the kinetic equation for the time evolution 
of the number density n^ (T) reads 



l"=<^)^54<»'-- 



~W:- 



3^12 



(18) 



The backward (inverse) direction is also valid. Note that n{T), Eq. ((T8)) . and n{p,T), Eq. (O, arc related with each 
other via n{T) = N{T)/V = g / {2t{^)^^ p{mi) j p^ dp niijj^T) and therefore n(p, T) is a Lorentz scalar whereas n{T) 
not. The thermal decay and production rate dW/dVdt have been discussed in j32| . In Boltzmann limit and assuming 
that the repulsive interaction does not contribute meaningfully to the overall relaxation time, the decay time in rest 
frame is given in textbooks. 
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Fig. 2: Left panel: ratio 77/^ vs. T. Right panel: thermal change of the relaxation time in Hagedorn fluid. The point-type and 
excluded-volume approaches are compared. Dashed-dotted line gives the speed of sound. 

where p is the pressure. / is a step functions for particle distinguishaility; 1 = 2 for indistinguishable and / = 1 for 
distinguishable particles. M is the hadronic reaction matrix. 

In right panel of Fig. [^ the relaxation time in fm, Eq. P^ . is given as function of T in GcV. We note that increasing 
the temperature T leads to reducing the relaxation time r. It might mean that the decay processes get faster when T 
increases. Near Tc, the effect of excluded-volume approach is considerable. 

As an application of these results, we mention the cosmological viscous models [33|, which require a complete set of 
thermo [l7| and hydrodynamic equations of state in order to solve the evolution equation in early universe and study 
the nucleosynthesis. 



IV. RATIO OF SHEAR VISCOSITY OVER ENTROPY DENSITY 



Using Anti de Sitter space/Conformal Field Theory (AdS/CFT) methods [34|. it has been argued that the ratio 
rj/s seems to have a universal lower bound in any physical system. The bound is rj/s > 1/Att. It has been found to 
this value saturated for a large class of strongly interacting systems with a dual description (string theory in anti-de 
Sitter space) [IJ, Is^- Also, at temperatures below Tc, i.e. Hagerdorn-type models, this ratio has been analyzed using 
chiral perturbation theory |36| , coupled Boltzmann equations of pions and nucleons in low baryon number densities 
Q , gas of massless pions [33, relativistic mean field models with scaled hadron masses and couplings [lal and viscous 
relativistic hydrodynamics [38|. In Fig. |31 wc draw ij/s calculated in the present model as a function of T. It is clear 
that ry/s starts from a much higher values than the AdS/CFT lower bound one. It comes closer to it with increasing 
T. 

Also, it is obvious that rj/s would be reduced with increasing the Hagedorn mass spectrum. It is a universal 
property that the quantities which are depending on thermodynamics are suppressed with increasing the mass. Shear 
viscosity 77 apparently follows this behavior, which can be realized primarily due to the enhancement of the massive 
resonances that leads to a decrease in the mean free path of a particle and a correspondin g iii crease in the average 
binary collision cross section. The results are compatible with the ones introduced in Ref. J15l | which are compared 
with the UrQMD simulations. Nevertheless, we notice however, that for Hagedorn resonance gas treated within the 
rate equation approach, an upper limit of rj/s is found to be as small as the KSS lower bound near Tr jlSll . 

The UrQMD model calculations give saturated 77/s-values in the region 0.125 < T > 0.166 GeV [I^l- Within this 
interval, 77 quantitatively equal to the entropy density s. Another difference between UrQMD and the present work 
is graphically shown in the T-region, 0.06 < T > 0.16 GeV, where it appears that the UrQMD model drastically 
underestimates the values of rj/s. It is clear that the massive resonances likely contribute dominantly, especially high 
T. The ratio 77/s gradually reduces. 

The effects of excluded volume on 77/s, which as discussed previously takes into consideration - at least - the van der 
Waals repulsive interactions, are also illustrated graphically in the same figure. The viscosity coefficients are enhanced 
due to this additional interaction. Finally, we notice that the values of ?//s very close to Tc are close to the upper 
bound of 0.24 obtained from viscous hydrodynamic calculations [331 of elliptic flow. 

To compare with the lattice QCD calculations below Tc [40|, we need to rearrange the configurations the Hagedorn 
fiuid to be compilable with the lattice QCD configurations. Since, it is not trivial to match the two configurations, we 
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Fig. 3: The shear viscosity to entropy density ratio 77/5 as a function of temperature T for the hadron resonance gas with 
Hagedorn mass spectrum. Horizontal dashed line give the AdS/CFT lower bound. 

think that a future study should be devoted to this subject. At higher T, the universal properties of bulk viscosity near 
the QCD phase transition are introduced in Ref. Q- On the other hand, there are several lattice QCD estimations for 
the lower bound [4l|. It is very essential to mention that these lattice QCD calculations have two limitations. First is 
the use of quench approximation, i.e. without quark pair creation-annihilation effects on vacuum, and second is the 
use of an ansatz for the spectral function. 
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